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Atomic oxygen adsorption on the Mo(112) surface has been investigated by means of first- 
principles total energy calculations. Among the variety of possible adsorption sites it was found 
that the bridge sites between two Mo atoms of the topmost row are favored for O adsorption at 
low and medium coverages. At about one monolayer coverage oxygen atoms prefer to adsorb in 
a quasi-threefold hollow sites coordinated by two first-layer Mo atoms and one second layer atom. 
The stability of a structural model for an oxygen-induced p{2 x 3) reconstruction of the missing-row 
type is examined. 

PACS numbers: 68.43.Bc, 68.43.Fg, 73.20.At, 73.30.-l-y 



I. INTRODUCTION 

Molybdenum surfaces attract substantial current in- 
terest because of their important applications and funda- 
mental properties. Molybdenum is one of the most inter- 
esting catalysts for its various oxidation states and a wide 
range of chemical reactivity. Oxidation of Mo strongly 
alters the dissociative adsorption of several molecules. 
Thus, oxygen adsorption on the low-index planes of bcc 
Mo has been widely studied using various experimental 
techniques. The preferred adsorption sites and ordered 
structures have been identified (see Ref. Q and references 
therein) . 

The Mo(112) surface although open, is also very sta- 
ble and is exploited as a substrate for the growth of or- 
dered silica layersSii or epitaxial growth of M0O2. The 
furrowed structure of Mo(112) has been used to study 
anisotropy of interactions in adsorbed submonolayers 
of alkali^ and alkaline-earth metal^ atoms. Recent ex- 
perimental findings, reporting that a presence of small 
amount of oxygen adatoms promotes a longer periodic- 
ity of low-coverage linear-chain-structures®'^ formed in 
the Mo(112) furrows, and responsible for commensurate 
phase transition in adsorbed Sr layers,* motivated this 
theoretical examination of O adsorption on Mo(112). 

According to experimental studies of low-coverage O 
adsorptionSiiSiii at Mo(112), oxygen atoms are expected 
to adsorb in trough sites between close-packed Mo rows 
and the sequence of most stable atomic oxygen sites can 
be ordered as: quasi-threefold hollow, long bridge, and 
on-top. 

Oxygen adsorption on the more open metal surfaces is 
known for long time to induce their reconstruction. One 
of the best known examples is the (2 x 1) added-row re- 
construction of the fee Cu(llO) surfaceiSiH which has a 
similar trough-and-row structure to the bcc (112). Dif- 
ferent possibilities for O and C induced reconstructions 
of Mo(112) and their effect on the surface band structure 
were discussed^ but it is only recently when Schrocder 
et a?ii*i^ reported results of a detailed experimental study 
of the transformation of the clean Mo(112) surface under 
oxygen exposure into an epitaxially grown M0O2 film. 



They found that an oxygen-induced p(2 x 3) reconstruc- 
tion precedes the M0O2 formation. This phase is very 
stable and passivates the Mo(112) against further oxida- 
tion in its initial stages. 

The principal goal of this investigation is the deter- 
mination of the oxygen modification of a Mo(112) sur- 
face and identification of preferred adsorption sites dur- 
ing low-coverage atomic oxygen adsorption, varying from 
a small fraction of a monolayer (ML) up to one mono- 
layer. The stability of different sites and configura- 
tions of O adatoms is explored including the observed 
p(l X 2) pairing- row overlayer structure which can be 
viewed as the one-dimensional M0O2 oxide structure. 
To the best of our knowledge such calculations have not 
been performed so far. We also investigate the stabil- 
ity of the structural model of the experimentally discov- 
ered oxygen induced (2 x 3) missing-row reconstruction 
of Mo(112),iii^ 

The remaining part of the paper is organized as follows. 
Section II outlines the methodology of calculations. It is 
followed by a presentation and discussion of our results in 
Sec. III. The paper is summarized in a concluding section. 



II. COMPUTATIONAL METHODS 

The calculations are based on density-functional 
theory ^^'^^ in the generalized gradient approximation — 
for the exchange-correlation energy functional. The 
electron-ion interaction is represented by the projector- 
augmented- wave potentials'^ as implemented^ in the Vi- 
enna ab initio simulation package (vASP)i^ The ground 
state is determined by solving the Kohn-Sham equations 
applying a plane waves basis set with energy cutoff of 
300 eV. A clean Mo(112) surface is modeled by periodic 
slabs consisting of seven molybdenum layers separated 
by empty space of sufficient thickness (eight equivalent 
Mo layers). Oxygen atoms are adsorbed in various sites 
on one side of the slab. The positions of Mo atoms in the 
four topmost layers, and all O atoms are fully relaxed un- 
til the forces on atoms converge to less than 0.02 eV/A. 
A dipole correction22i2^ was applied to compensate for 
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asymmetry of the potential at the two sides of the slab. 
To sample the Brillouin zone (BZ), Monkhorst-Pack It- 
point meshesS^ were applied, corresponding to 20 x 12 x 1 
k-points in the 1x1 primitive unit cell, and accordingly 
reduced for larger cells. To improve the quality of BZ 
integrations for the fractional occupancies a Methfessel- 
Paxton method-^ with a broadening of 0.2 eV is applied. 
In the calculations for slabs the spin-polarization effects 
were not included. The binding energy of oxygen is cal- 
culated with respect to one-half of the energy of a spin- 
polarized O2 molecule placed in a large rectangular unit 
cell. 

Our slab reproduces well^ both the electronic proper- 
ties (the work function equals to 4.12 eV) and the ge- 
ometric structure of the upper layers of bare Mo(112), 
including a 17% contraction of the interplanar distance 
for the topmost atomic plane, and a relatively small ver- 
tical relaxation of the second (—0.9%), third (3.0%), and 
fourth (—0.3%) interlayer spacing. 



III. RESULTS AND DISCUSSION 
A. O adsorption on unreconstructed Mo(112) 

In order to determine the preferred sites for O adsorp- 
tion on a clean Mo(112) we consider first very low cov- 
erages of 1/8 ML corresponding to a single O atom in 
a 4 X 2 surface unit cell. We define the coverage Q as 
the ratio of the number of adatoms to the number of 
substrate atoms in a surface unit cell. Following exper- 
imental findingsiSiii we explore several possible adsorp- 
tion sites for O atoms: long-bridge, atop, quasi-threefold 
hollow, and short-bridge. The locations of these sites are 
shown schematically in Fig. 

At this coverage the distances between adatoms are 
large (10 A across, and 11 A along the rows), and the 
0-0 interaction is negligible. The strongest binding for 
this coverage is found when the O atoms are adsorbed 
in short bridge sites in a topmost Mo row (Table [}. 
Note that the binding (adsorption) energy, which is calcu- 
lated with respect to the energy of one-half of the energy 
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FIG. 1: Top view of the upper layers of the ideal furrowed 
Mo(112) surface and the possible locations of adsorption sites 
for O atom. From left to the right two possibilities for quasi 
threefold hollow (hoi), short bridge (bri), long bridge (lb), and 
atop (top) sites are represented by black circles. 



TABLE L Binding energies of O atom adsorbed in short 
bridge positions at a topmost Mo rows and quasi-threefold 
hollows of Mo(112) in different adsorbate structures. The last 
line entry corresponds to the bottom-right (1x2) pairing- row 
structure displayed in Fig. |5| 
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of a free O2 molecule, from the total energy difference 
of the systems with and without adsorbate, is defined 
here as positive»i A second stable site at low coverages 
is the quasi-threefold hollow site situated between two 
row atoms and one trough atom with the binding en- 
ergy about 120 meV lower. In the following, for brevity, 
the former sites are referred to as bridge and the lat- 
ter as hollow sites. The other sites appear either un- 
stable (long-bridge in the trough, and atop Mo atom in 
the rows — the O atoms shift from there to hollow or 
to short-bridge sites, respectively) or show a distinctly 
weaker binding (atop is 0.6 eV weaker bound) compared 
to the strongest bound O in short-bridge sites. Interest- 
ingly, the in-trough located short-bridge sites are equally 
unfavored as on-top sites. Actually, they both end up in 
a shifted short-bridge site which could be called a hollow 
site of second type. This site is coordinated by one first- 
layer atom and two second layer Mo atoms. A preference 
for the bridge sites on the top Mo row is somewhat un- 
expected because oxygen atoms usually tend to occupy 
a higher coordinated sites on transition metal surfaces. 
Therefore, it contradicts the experimental models which 
suggest a quasi-threefold hollows as most favored sites for 
low-coverage 0-adsorption on Mo and W(112)Ai2iS& It 
also differs from O adsorption on the furrowed Cu(llO) 
surface where both experiment and theorjii^^ predict high- 
coordinated hollow sites as most favored. 

Table in presents the calculated binding energies. The 
O coverage was increased either by placing additional O 
atoms in the surface unit cell or by reducing cell size. 
Test calculations performed at the coverage Q = 0.5 for 
the bridge- and hollow-site p{l x 2)-0 structures in 1 x 2 
and 2x2 unit cells, show that the binding energies agree 
within 15 meV and 13 meV, respectively. Hence, the 
maximum error bar in the determination of the relative 
energies, with respect to different size of unit cells and 
different number of k-points applied, is estimated as ±15 
meV. For a given coverage, we compare the energies of 
different patterns calculated in a unit cell of same size 
and the error is much below this limit. The binding en- 
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FIG. 2: Schematic view of the adatom structures formed at 
various O coverages on the Mo(112) surface. The configu- 
rations of the upper row correspond to oxygen adsorbed in 
bridge sites while those of the lower row to O adatoms in hol- 
low sites. For all the displayed O configurations both types of 
adsorption sites are considered (TablePi except of the bottom- 
right pairing-row structure which appears only for the hollow 
sites. 



ergy differences for low-coverage structures amount to at 
least several tens of meV, thus they are well beyond the 
above error bars and different structures can be easily 
discriminated. 

The results of TablelHseem to agree with the scenario of 
an initial stage of O adsorption based on the low energy 
electron diffraction (LEED) experiment.^ According to 
that, beginning from lowest coverages, the oxygen atoms 
tend to form atomic rows along the [111] direction with 
twice the spacing between atoms in the close-packed Mo 
rows. The distance and phase of adjacent O rows are 
initially random. At increased coverage 0-rows tend to 
occupy in-phase every second row. This results in for- 
mation of a (2 X 2) pattern at 0.25 ML (Fig. |2J|. At 0.5 
ML coverage, the sites along each Mo row are occupied 
and the p{2 x l)-0 structure is developed. This pattern 
is more stable than the c(2 x 2) which is characterized by 
antiphase 0-rows or the p(l x 2) structure with O atoms 
in each site along every second row. It is also clear that 
for coverages 9 < 0.5 ML the bridge sites remain most 
favored. 

With increased coverage (9 = 0.5), when the O atoms 
begin to occupy neighboring sites along the rows [p{l x 2) 
structure] , quasi-threefold hollow sites become more pre- 
ferred than bridge sites (Table Q}. Note, however, that 
for this coverage the p{2 x 1) and c(2 x 2) structures 
with a more homogeneous distribution of O atoms still 
remain energetically most favored. These results suggest 
that 0-0 repulsion is stronger for bridge than for hollow 
adsorption. In the latter positions the O adatoms are 
immersed in the electron charge distribution of the Mo 
substrate which screens the dipole-dipole repulsion. For 
a higher coverage our results predict a structural phase 
transition, from a p{l x 1) phase with adatoms preferen- 



TABLE II: The geometry of oxygen adsorption on Mo(112). 
Adi2 is the change in the vertical relaxations of the topmost 
interlayer distance, /ii{2) is the height of O atom with respect 
to the topmost (second) Mo layer, / is the O-Mo bond length 
(distance). For a clean Mo(112) Adi2 = -17.0%. 
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tially occupying the bridge sites, to a p(l x 1) structure 
where adatoms occupy hollow sites. The energetic differ- 
ences for O adsorbed in hollow sites at different coverages 
and structures are smaller than for O in bridge-top sites. 

For both types of energetically favorable sites a de- 
crease in the binding energy with increasing coverage is 
clearly seen (Tabled). While for hollow sites the changes 
in the binding energy do not exceed 190 meV, for the 
whole range of coverages considered they are four times 
larger for bridge sites. In the latter case, they increase to 
0.7 eV for 9 > 0.5 in the cases when the nearest neigh- 
bor sites parallel to the troughs are occupied. For bridge 
sites the stability of their structures enhances with align- 
ing adatoms normal to the furrows whereas for hollow 
sites the opposite can be observed. 

Our results show (Table P) that the (1 x 2)-0 pairmg- 
row reconstruction, illustrated in the bottom-right panel 
in Fig. 121 where O atoms occupy hollow sites on both 
sides of the Mo row is equally favored as the (1x1) struc- 
ture with O atoms homogeneously distributed in hollow 
sites on one side of the Mo rows. In LEED experiments 
one originally observes (1x1) patterns, and a (1x2) pat- 
tern is formed only after annealing to ca. 500-600 Ki2ii£ 
The calculated barrier (0.63 eV) for O diffusion from a 
hollow site on one side of the trough to the one located 
on the other side (Fig. ^ confirms that this configuration 
can be reached only during activated adsorption process. 
The stability of this structure means that atoms of a Mo 
row screen very effectively the repulsive interaction be- 
tween O atoms of the pairing rows. 

The changes in the average geometry of Mo(112) lay- 
ers with O coverage are presented in Table ITIl It is seen 
that O adatoms very differently influence the geometry 
of Mo(112). The O atoms adsorbed in bridge sites have 
a stabilizing effect on the Mo(112) surface which is man- 
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FIG. 3: Work function changes induced by O atoms adsorbed 
on the Mo(112) surface in short bridge positions (circles) of 
the first layer and threefold hollows (diamonds). The work 
function for a pairing-row p(l x 2)-0 structure is marked by a 
square. Triangles down and up (at 0.5 ML) label the entries 
for the c(2 x 2) and (1 x 2) structures, respectively. 



ifested in a diminished contraction (derelaxation) of the 
interlayer distance between the topmost Mo layers with 
increasing O coverage. In contrast, oxygen adsorbed in 
the quasi-threefold hollow sites enhances relaxation of the 
topmost layers. For higher coverages the vertical shifts of 
the atoms of the particular planes relative to the center- 
of-mass position (surface rumpling) are much smaller or 
negligible. Interestingly, the pairing-row p(l x 2)-0 struc- 
ture involves some parallel shifts of top Mo rows. Every 
second Mo row decorated by O atoms from both sides, 
shifts in the [TTl] direction by about 0.3 A, whereas the 
other, "bare" Mo-row, not coordinated directly by oxy- 
gen, shifts in an opposite direction by approximately 0.1 
A. Furthermore, the bare Mo-rows are shifted down (by 
~ 0.09 A), while the other with coordinating 0-rows, 
shift up by the same amount with respect to the center- 
of-mass position. A single O atom adsorbed in a bridge 
site in a 4 X 2 cell induces substantial vertical shifts of 
the two Mo atoms closest to oxygen (0.1 A up). The re- 
maining atoms of the same Mo row do not change their 
positions. For an O atom adsorbed in a quasi-threefold 
hollow, one of the two coordinating Mo atoms of the up- 
per row moves down (0.1 A) while the other shifts 0.1 A 
upwards, and the O atom shifts laterally (by 0.08 A) in 
the direction of the latter. All these adjustments enhance 
the coordination of the oxygen atom. No change in the 
vertical alignment of the Mo rows is observed. The de- 
crease in the binding energy with coverage is, in general, 
accompanied by a small, though not regular, decrease in 
the vertical distance of the O adatom to the topmost Mo 
plane. 

The distance /12 to the second Mo layer increases with 
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FIG. 4: Atom resolved density of states for O atom and the 
nearest Mo atoms of the topmost two layers of Mo(112) at 0.25 
ML of oxygen coverage. The left column displays the states 
for O adsorbed in bridge site while the right one in hollow 
site. The energy is measured with respect to the Fermi level. 



coverage for bridge sites and decreases for hollow sites. 
For the bond length I an opposite trend is observed. Vari- 
ations of I correlate with a reduced contraction of the 
topmost Mo-layers spacing. 

The calculated work function variation due to adsorp- 
tion of electronegative O exhibits a strong increase with 
coverage which is consistent with a general picture of 
oxygen chemisorption on transition metals. Note that 
this increase is almost three times larger for O atoms 
adsorbed in short-bridge positions than for those in the 
hollow sites. This is presumably due to a much stronger 
polarization of protruding O adatoms in bridge than hol- 
low sites. The work function variation seems to be con- 
sistent with experiment which reports a rapid increase 
of work function followed by a slower decrease for higher 
oxygen exposures. Our binding energy data (Table P) 
show that for 0.5 < coverages < 1 ML a site exchange 
transition from bridge to hollow sites occurs, therefore a 
work function decrease marked schematically in Fig.|2|by 
a dashed line should accompany this transition. 

Figures 01 and |S1 display the local partial density of 
electronic states (LDOS) at the O atom and the neigh- 
boring atoms of the two topmost layers, for a low and a 
full monolayer coverage with O atoms in bridge and hol- 
low sites. The main contribution to the bonding comes 
from the O 2p orbitals. For 1/4 ML coverage, with O in 
bridge sites this contribution consists of two partly over- 
lapping peaks, with the higher of them located at —4.5 
eV. For O in hollow site the region of 2p states bonding 
is narrower (~ 1 eV wide) with only one peak centered 
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FIG. 5: Same as in Fig.^Jfor 1 ML coverage. 

at —5.5 eV. As is seen from Fig. ^ the O 2p orbitals 
hybridize strongly with the Mo Ad states. For oxygen 
adsorbed in the bridge site the bonding is hmited only to 
the Mo atom of the topmost row while for O adsorbed 
in hollow sites both the Mo atom of the upper row and 
the one of the trough contribute to the bonding. The 
density of states of the Mo atom of the trough resembles 
very much that of Mo bulk^^ which is dominated by a 
conduction band formed of localized 4c? orbitals. 

For a 1 ML coverage (Fig. El, the regions of bonding 
widens both for bridge and hollow sites adsorption. It 
extends from —5.8 to —3.5 eV with two peaks at about 
—5.5 and —4.2 eV, for the bridge site adsorption. Again, 
the O 2p orbitals hybridization with the Mo id orbitals 
contributes mostly to the bonding. For O adsorption in 
hollow sites, the single peak due to O 2p orbitals is split 
into three distinct components and the range of bonding 
has a similar width to that for adsorption in bridge sites. 
The hybridization of the 5s states of the upper Mo row 
with O 2p orbitals is also increased. 

B. Missing-row reconstruction of Mo(112) 

It has been recently reported that O adsorption in- 
duces a p{2 X 3) surface reconstruction of the missing 
row typeii^ A structure model of reconstruction was 
proposed [Fig. El^a)], resulting from a combined scan- 
ning tunneling microscopy and LEED study. In order 
to check stability of this structure, which assumes oxy- 
gen placed both in hollow and on-top Mo-atom positions, 
we have used it as an input model for the total energy 
calculations. The structure optimization leads to the O- 
pattern [Fig. El^b)] which differs from the one suggested 
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FIG. 6: A p(2 X 3) unit cell of the missing row reconstruction 
of the Mo(112) surface induced by O adsorption, (a) experi- 
mental model, (b) structure optimized positions of model (a), 
an example of the calculated more stable structure (c) . Oxy- 
gen atoms are represented by small dark spheres. 

by experiment In the optimized structure, the oxygen 
atom placed originally atop the Mo atom of the second 
Mo layer is shifted to the quasi-threefold hollow, thus 
aligning the O atoms with a Mo row. This is not too 
surprising in view of our findings, discussed in Sec. IIIA, 
that on-top sites belong to the least favored ones for O 
adsorption at Mo(112) and are 0.6 eV less stable than 
hollow sites. It means that an O atom sitting in a hol- 
low site provides a more effective coordination to the Mo 
atoms of second layer than atop configuration. 

The calculated p{2 x 3) structure belongs to the most 
stable ones among several considered configurations of O 
atoms in the unit cell, but even limited search for the 
total energy minimum, within a class of different O as- 
signments, allows to find other, more stable structures, 
involving oxygen atoms of only hollow sites. An example 
of such a structure (prefered by 446 meV per unit cell) is 
given in Fig. |Sfc). Thus, though the optimized structure 
of Fig.EIb) corresponds best to the one suggested by the 
experimental studies^'"'^^ it does not provide a global en- 
ergy minimum and should be considered as one of the 
most likely but metastable configurations. 



IV. CONCLUSIONS 

The systematic total-energy calculations of atomic 
oxygen adsorption on the relaxed Mo(112) surface are 
presented and discussed. In contrast to experimental as- 
signments our calculations identify the short-bridge sites 
of the uppermost Mo row as the most stable for O ad- 
sorption at lower coverages. They are by about 120-150 
meV per atom more stable than the quasi-threefold hol- 
lows. The binding energy of O atoms adsorbed in the 
stable bridge and quasi-threefold hollow sites decreases 
with increasing coverage. Our results show that O atoms 
very differently modify the spacing between the topmost 
two Mo(112) layers. Thus, measurements of lattice re- 
laxation can provide a method for identification of the 
adsorbed oxygen sites. At one monolayer coverage a site 



6 



exchange transition from bridge to quasi-threefold hoUow 
sites occurs. The calculated work function variation with 
O adatom coverage is found to be consistent with recent 
measurements. The stability of the p{l x 2) pairing-row 
reconstructed O pattern predicted by LEED experiment 
is confirmed. The p{2 x 3) missing-row reconstruction is 
investigated. The stability of the configuration suggested 
by experiment requires slightly modified positions of O 
atoms. 
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